Airy beams from a microchip laser 
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It is theoretically shown that an end-pumped microchip laser formed by a thin laser crystal with plane-plane 
but slightly tilted facets can emit, under appropriate pumping conditions and near a crystal edge, a truncated 
self-accelerating Airy output beam. © 2013 Optical Society of America 
OCIS codes: 140.3410, 350.5500, 190.4420, 140.3580 
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Microchip lasers are an important class of monolithic 
solid-state lasers which have attracted a considerable 
interest as single-frequency and ultracompact devices 
[1-3]. A microchip laser generally consists of a thin rare- 
earth-doped laser crystal with dielectric coatings on both 
faces to form a plane parallel cavity, end-pumped by a 
diode laser. Owing to the planar nature of the cavity, the 
transverse modes of a microchip laser are not determined 
by mirror curvatures, rather by weaker, residual guiding 
effects, such as gain or thermally-induced guiding [4-9]. 
In this Letter we consider a microchip laser with two pla- 
nar but slightly tilted facets, and show that end-pumping 
near the edge of the laser crystal can lead to laser emis- 
sion into self-accelerating Airy beams, a class of opti- 
cal beams which has received a great interest in the re- 
cent few years [10-17]. Airy beams show rather unique 
properties, such as diffraction free propagation along a 
curved (parabolic) path and self healing, i.e., restoring 
their canonical form after passing small obstacles. Airy 
beams found recently important applications in optics, 
including optical trapping [12], plasma waveguiding [15] 
and nonlinear frequency generation [16,17]. 
Let us consider a microchip laser operating at the wave- 
length A made of a thin crystal (length refractive in- 
dex Us) with two highly-reflecting and dielectric-coated 
planar facets [Fig. 1(a)]. For the sake of simplicity, we 
first consider the case of a one transverse spatial coor- 
dinate. The two crystal facets are not perfectly paral- 
lel each other, but are tilted to form a small angle 
The microchip is end-pumped by a Gaussian beam, fo- 
cused at a distance px from the upper edge of the crys- 
tal [see Fig. 1(a)]. Field propagation in the cavity can 
be described within an effective-index and mean-field 
model [9], in which the intracavity field envelope is 
assumed to slightly vary in one cavity round trip. Prop- 
agation back and forth between the two mirrors is equiv- 
alent to propagation along a lensguide with a complex 
refractive index n{x) = + An(x) n^, where the 
real and imaginary parts of An(x) account for (weak) 
index and gain guiding effects, respectively [9]. Let V^(x) 
be the field envelope at the reference plane 2: = in 
the cavity; after a cavity round-trip, the field envelope is 
transformed into the one ip'{x) = Tip(x)^ where T is the 



round-trip propagation operator given by 

T = exp(— 7) ex.p{—iklH) ex.p{—iknsaxx) ex.-p{—iklH). 

(1) 

In Eq.(l), /c = 27r/A is the wave number, 7 is the single- 
pass logarithmic loss of the cavity [18], and 
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Note that tlj\x) is obtained from ?/^(x) by application 
of the four operators Ti = exp(— 7), T2 = exp(— i/c/ff), 
T3 = ex.p{—iknsaxx) and T4 = T2 describing the effects 
of cavity loss, forward propagation in the crystal, phase 
shift introduced by the mirror tilt, and backward propa- 
gation in the crystal, respectively. In the mean- field limit, 
one can assume T 1 — 7 — 2iklH — ikngO^xX. Indicat- 
ing by Tr = 2lns/co the round-trip transit time in the 
cavity, where cq is the speed of light in vacuum, the tem- 
poral evolution of the field envelope ?/^(x, t) is then given 
by dt^{x,t) [ij\x) -iIj{x)]/Tr, i.e. 
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where Tc = Tr/(2^) is the cavity photon lifetime. Ne- 
glecting thermal lensing effects [7] and assuming that 
the laser operates at the resonance frequency of the gain 
medium, so that gain-related index guiding effects can 
be neglected as well [7,9], the refractive index An in 
Eq.(3) can be taken to be purely imaginary and given by 
An{x) = —i(jgN{x)/(2k)^ where (jg is the cross-section 
of the lasing transition and N{x) is the (saturated) pop- 
ulation inversion distribution in the crystal. For a Gaus- 
sian pump beam with a Rayleigh range larger than the 
crystal thickness / and for the laser operated near thresh- 
old, for which saturation effects can be neglected [7], one 
has N{x) — No{x) = Np ex.p[—2{x — px)'^ /w^]^ where Wp 
is the pump beam spot size. Above threshold, in the 
mean-field limit and neglecting spatial hole burning, fol- 
lowing Ref. [8] saturation effects can be included in the 
model by taking N{x) = A^o(^)/(l + IV^P), where the in- 
tracavity field intensity is normalized to the saturation 
intensity of the lasing transition. The on-axis single-pass 
(unsaturated) photon gain coefficient Qp in the medium 
is related to the (unsaturated) peak population inversion 
Np by the simple relation Qp = (JgNpl [18]. Crystal trun- 
cation at X = introduces an effective barrier potential 
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Fig. 1. (Color online) (a) Schematic of an end-pumped 
microchip laser with tilted mirror, (b)-(d) Numerically- 
computed transverse mode profile at lasing threshold 
(solid curve) for decreasing values of the logarithmic 
cavity loss, showing the transition from an asymmetric 
Gaussian beam to a truncated Airy beam: (b) 7 = 0.142, 
(c) 7 = 0.0651, and (d) 7 = 0.019. The gain threshold 
gp = 9pth is Qpth = 0.326 in (b), Qpth = 0.2172 in (c), 
and Qpth = 0.1086 in (d). The dotted curve in the figures 
is the profile of the Gaussian population inversion ATVq, 
whereas the dashed curve in (d) is the intensity profile 
of the Airy mode |V^nP with = 23.224. (e) Free-space 
propagation of the output laser beam of Fig. 1(d) [inten- 
sity distribution in the (x, z) plane], (f) Same as (d), but 
for a gain parameter Qp two times (left panel) and three 
times (right panel) above threshold. 



for the paraxial lasing field because of total internal re- 
flection, which is modeled by requiring = for x < in 
Eq.(3). As discussed below, such an edge effect is of ma- 
jor importance for the laser to oscillate on a (truncated) 
Airy transverse mode. Introducing the normalized spa- 
tial and temporal variables = x/L and r = t/T^ where 

L=(^-\'\ T = ^('^\'\ (4) 

yn'laxk'^ J ' Co \alk J 



Eq.(3) can be cast in the dimensionless form 

where F = F{£^) = exp[— 2(^L — px)"^ /ujp] is the normal- 
ized profile of the unsaturated population inversion in 
the crystal and gp = (JgNpl the peak gain coefficient. 
The thresholds and profiles of the various transverse 
modes sustained by the microchip laser can be computed 
by numerical analysis of eigenvalues and eigenfunctions 
of the (linear) operator T by taking AN{x) = ANo{x). 
In particular, the lowest-order lasing mode can be 
computed by a standard Fox-Li iterative method [19]. 
The profile of the lasing mode turns out to depend 
mainly on the two parameters {T/Tr)^ and px/L^ 
i.e. by the logarithmic cavity loss, tilting angle, and 
distance of the pump beam from the crystal edge. This 
is because two different physical mechanisms compete 
in the formation of the transverse laser mode. On the 
one hand, gain guiding provided by the pump tends to 
localize the lasing mode near the gain region, as in usual 
microchip lasers with parallel facets [4-8]; on the other 
hand, the mirror tilting pushes the light beam toward 
the edge of the crystal, where it is back reflected owing 
to total internal reflection, thus providing an effective 
trapping of light which was discussed in [20] in another 
optical setting. Depending on the strength of these two 
competing processes, two distinct transverse confining 
regimes can be found. When p/L^ is very large (i.e. 
the crystal is pumped far from the crystal edge) and 
the cavity losses are not too small (i.e. the gain guide 
is strong enough), the gain-guiding mechanism prevails 
and the lasing mode looks like a slightly asymmetric 
Gaussian beam, similar to the off-axis Gaussian modes 
with complex arguments studied in [21,22]. Conversely, 
if the crystal is pumped closer to the edge of the crystal 
and the logarithmic cavity loss 7 is sufficiently small, 
the edge trapping mechanism prevails. In this case, gain 
and loss terms entering on the right hand side of Eq.(5) 
can be treated as perturbation terms, and the transverse 
mode profiles of the cold cavity are defined by the Airy 
equation — (^i/j = ai/j with the boundary condition 
V^(0) = 0, i.e. V^n(0 = ^i(^ + ^n), where is the nth 
zero of the Airy function and <j = is the eigenvalue 
that determines the transverse mode frequency. Hence, 
the lasing mode will be determined by the (truncated) 
Airy eigenmode V^n(0 ^i^^ largest overlap with 
the pump beam profile. As an example, Figs.l(b-d) 
show the transition of the lasing mode at threshold 
from an asymmetric Gaussian-like beam to a truncated 
Airy beam. In the simulations, different crystals were 
considered, with decreasing values of cavity loss 7 
(controlled e.g. by the mirror reflectance) and equal 
distance px of the pump beam from the crystal edge. 
The lowest-order (lasing) mode at threshold has been 
computed by the Fox-Li iterative method using the 
exact round-trip propagator T, defined by Eq.(l). Pa- 
rameter values used in the numerical simulations refer to 
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a Nd:YV04 microchip laser [7,8] (n^ = 2.17, A = 1064 
nm) with crystal thickness / = 400 /im and tilt angle 
ax = 0.5 mrad, which correspond to spatial and tempo- 
ral scales L \1 fim and T 53 ps [see Eq.(4)]. The 
pump distance from the crystal edge is px = 373 /im, 
whereas the pump spot size is Wp = 17 /im. As for 
large cavity losses the gain guiding is the dominant 
mechanism of transverse confinement and the lasing 
mode looks like a slightly asymmetric Gaussian beam 
[Fig. 1(b)], at lower cavity losses index guiding provided 
by edge refiection and mirror tilt prevails and the lasing 
mode is well approximated by a truncated Airy beam 
with its largest lobe overlapped with the pump beam 
[see Fig. 1(d)]. The free-space propagation of the output 
beam in this case is depicted in Fig. 1(e), clearly showing 
the characteristic parabolic (accelerated) path of the 
Airy beams. It should be noted that the transition from 
Gaussian-like to Airy-like modes shown in Fig.l could 
be observed by assuming different microchip lasers with 
equal loss 7 but with increasing tilt angle Laser 
emission in the Airy-like mode persists for the microchip 
operated well above threshold, where saturation effects 
can not be neglected. As an example. Fig. 1(f) shows the 
numerically-computed steady-state lasing modes, which 
account for saturation of the population inversion, 
for the same conditions of Fig. 1(d) but for the laser 
operated with a gain parameter gp two and three times 
above threshold. As one can see, gain saturation does 
not appreciably modify the transverse mode profile. 
The previous analysis can be extended to the case of 
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Fig. 2. (Color online) Two-dimensional transverse in- 
tensity profile \ilj{x^y)\'^ of the lasing mode at threshold 
for (a) 7 = 0.0632, (b) 7 = 0.0298, and (c) 7 = 0.0106 
(other parameter values are given in the text). The gain 
threshold is gpth = 0.354 in (a), gpth = 0.236 in (b), and 
gpth = 0.118 in (c). 

two- trans verse spatial variables x and y. Indicating by 
ax and ay the mirror tilting angles with respect to the 
X and y axes, in terms of the dimensionless variables 
^ = x/Lx^ T] = y/L and r = t/T, where L and T are 
defined by Eq.(4), the evolution equation of the field en- 
velope in the mean-field limit reads drip = Qip, where 
Q = i[dl+d'^-^-{ay/a,)rj] + iT/Tn)[gpF/il + \i,\^)-j], 
and F = F(^, r]) is the normalized profile of the 
population inversion in the crystal. For a circular 
Gaussian pump beam with spot size Wp focused at 
the point (px^Py) in the transverse plane, one has 
Fi^,rj) = exp[-2(ei: - p.f/wl - 2{r,L - Pyf/wl]. 



Provided that the gain guiding is weaker than the index 
guiding realized by the mirror tilt and refiection at the 
edges, the output laser mode is expected to be given by 
a two-dimensional truncated Airy beam, with its main 
lobe mostly overlapped with the Gaussian pump beam. 
As an example. Fig. 2 shows the numerically-computed 
two-dimensional distributions of the most unstable laser 
mode in a Nd:YV04 microchip laser, with a transition 
from a Gaussian-like to a two-dimensional truncated 
Airy profiles, for decreasing values of the logarithmic 
cavity losses 7 and for / = 400 /im, ax = ocy =0.2 mrad, 
Px = Py = 184 /im, Wp = 23 /im. 

In conclusion, it has been shown that microchip lasers 
with slightly tilted facets can emit transverse output 
beams in the form of truncated Airy beams. Such a 
result could be of interest for a direct generation of Airy 
beams from an ultracompact microchip device. 
Work supported by the Italian MIUR (Grant No. 
PRIN-2008-YCAAK). 
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